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ABSTRACT 

A model for the infrared emission of the high redshift uhraluminous in- 
frared galaxy FSC 10214+4724 is presented. The model assumes three com- 
ponents of emission: a dusty torus viewed edge-on, clouds that are associated 
with the narrow- line region and a highly obscured starburst. It is demonstrated 
that the presence of clouds in the narrow-line region, with a covering factor of 
17%, can explain why the mid-infrared spectrum of FSC 10214-f 4724 shows a 
silicate feature in emission despite the fact that its torus is viewed edge-on. It 
is also shown that the same model, but with the torus viewed face-on, predicts 
a spectrum with silicate emission features that is characteristic of the spectra 
of quasars recently observed with Spitzer. 

Key words: galaxies: active - galaxies: individual (FSC 10214+4724) - 
infrared: galaxies - dust: - radiative transfer: 



1 INTRODUCTION 

The z=2.285 galaxy FSC10214+4724 was discovered 15 
years ago by Rowan- Robinson et al. (1991). At the time 
of its discovery it was thought to be one of the most 
luminous objects in the Universe. It was later found 
to be gravitationally lensed by a factor 10-100 (Broad- 
hurst & Lehar 1995, Graham & Liu 1995, Serjeant et al. 
1995, Eisenhardt et al. 1996). It is still the best stud- 
ied ultraluminous infrared galaxy at that redshift. Rest 
frame optical-UV spectroscopy and spectropolarimetry 
showed that FSC 10214+4724 exhibits similar char- 
acteristics to Seyfert 2 galaxies (Elston et al. 1994, 
Lawrence et al. 1993, Goodrich et al. 1996). A large 
quantity of molecular gas is suggested by CO observa- 
tions (Brown & Vanden Bout 1991, Solomon et al. 1992, 
Scoville et al. 1995). The object is also bright in the sub- 
millimetre. It is therefore likely that FSC 10214+4724 
is similar to some local ultraluminous infrared galaxies 
which harbour simultaneously a luminous starburst and 
an AGN (e.g. Genzel et al. 1998). 

Radiative transfer models of emission from dust 
in a starburst or an AGN were used in an attempt 
to understand the origin of the infrared emission of 
FSC10214+4724. Rowan-Robinson et al. (1993) showed 
that good fits to the overall spectral energy distribu- 
tion (SED) could be obtained with both a starburst 



and an AGN model. The dust density in the latter 
falls off as and varies exponentially with azimuthal 
angle (Efstathiou & Rowan- Robinson 1995). Green & 
Rowan- Robinson (1996) concluded that starburst or 
AGN activity alone can not explain the SED of FSC 
10214+4724 but instead a combination of the two is 
needed. Granato, Danese & Franceschini (1996) found 
that a good fit to the SED of FSC 10214+4724 and other 
hyperluminous infrared galaxies could be obtained with 
pure dust-enshrouded AGN models. The limited wave- 
length coverage in the infrared is the main reason why 
modeling attempts in the mid-90s remained inconclu- 
sive. 

Recently Teplitz et al. (2006) presented a mid- 
infrared spectrum for FSC 10214+4724 which was ob- 
tained with the Infrared Spectrograph onboard Spitzer 
(Houck et al. 2004). This spectrum places strong con- 
straints on models of the infrared emission of this ob- 
ject. What is interesting about the spectrum presented 
by Teplitz et al. is that the silicate feature appears in 
emission. This is surprising given the type 2 nature of 
the AGN in FSC 10214+4724. Local examples of type 
2 AGN (e.g. NGC1068) show silicate features in ab- 
sorption in agreement with radiative transfer models of 
dusty tori viewed almost edge-on (Pier & Krolik 1992, 
Granato & Danese 1994, Efstathiou & Rowan- Robinson 
1995). The spetrum presented by Teplitz et al. shows no 
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signs of PAH features except possibly a weak 6.2^m fea- 
ture. This implies that the starburst if present is weak 
or highly obscured. 

In this letter radiative transfer models of dusty 
AGN and dusty starbursts are used to fit the now de- 
tailed SED of FSC 10214-^4724. Section 2 describes 
briefly the radiative transfer models to be used. In 
section 3 the results are presented and briefly dis- 
cussed. A flat Universe is assumed with A = 0.73 and 
Ho=71Km/s/Mpc. 



2 RADIATIVE TRANSFER MODELS 

2.1 Starburst models 

Radiative transfer models for starbursts were pre- 
sented by Rowan- Robinson & Crawford (1989), Rowan- 
Robinson & Efstathiou (1993), Kriigel & Siebenmor- 
gen (1994), Silva et al. (1998), Efstathiou et al. (2000), 
Siebenmorgen et al (2001), Tagaki et al. (2003), Do- 
pita et al. (2005). For this study the models of Efs- 
tathiou et al. (2000) which were extended to higher op- 
tical depths by Efstathiou & Siebenmorgen (2005) will 
be used. The models incorporate the stellar population 
synthesis models of Bruzual & Chariot (1993), a sim- 
ple model for the evolution of the molecular clouds that 
constitute the starburst and detailed radiative trans- 
fer that includes the effect of transiently heated par- 
ticles/PAHs (Siebenmorgen & Kriigel 1992). In Fig. 1 
the model predictions for starbursts of different opti- 
cal depths are compared. As expected the ratio of far- 
to mid-infrared luminosity increases with optical depth 
whereas the submillimeter spectrum remains largely un- 
affected. The model with ry ~ 50 predicts an SED sim- 
ilar to that of the nearby starburst galajcy M82 whereas 
the most optically thick model gives an SED similar 
to that of the heavily obscured ultraluminous infrared 
galaxy Arp220. 

2.2 AGN torus models 

Models of the infrared emission from dusty tori and discs 
were presented by Pier & Krolik (1992, 1993), Granato 
& Danese (1994), Efstathiou & Rowan-Robinson (1995), 
van Bemmel & DuUemond (2003), DuUemond & van Be- 
mmel (2005). Discussion centred on the problem of mini- 
mizing the silicate emission features at lOfim which were 
not observed in the ground-based 8-13^m observations 
of predominantly Seyfert galaxies by Roche et al. (1991). 
The non-detection of silicates was attributed either to 
the torus geometry (Pier & Krolik 1992, Efstathiou 
& Rowan-Robinson 1995) or to the destruction of sil- 
icate dust in the inner torus (Granato & Danese 1994). 
Rowan- Robinson (1995) and more recently Nenkova et 
al. (2002) proposed that the non-detection of silicates 
is due to the dumpiness of the dust but DuUemond & 
van Bemmel (2005) showed that dumpiness does not 
necessarily eliminate the silicate features. 
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Figure 1. Spectral energy distributions of Efstathiou et al. 
starbursts in which the star formation rate decays exponen- 
tially with time (time constant SOMyrs). The SED is shown 
for three different values of the initial optical depth in the V 
band of the molecular clouds that constitute the starburst. 
The initial star formation rate is assumed to be 1 Mq/ht 
and the starburst age is assumed to be 40Myrs. 



Recent results from Spitzer (Siebenmorgen et al. 
2005, Hao et al. 2005) showed that quasars emit emis- 
sion features both at 10 and 18fj,m. Although a detec- 
tion of silicate emission in a LINER was also reported 
by Sturm et al. (2005) it seems that silicate emission is 
found mainly in high luminosity AGN. This is an ob- 
servational fact that models of the infrared emission of 
AGN must therefore explain. The models should also 
explain the fact that some luminous type 2 AGN like 
FSC 10214-1-4724 (see also Sturm et al 2006) also show 
emission features at 10/^m. 

Another observational fact is that about half of 
the mid-infrared emission of the prototypical Seyfert 
2 galaxy NGC1068 is extended and coincides with the 
narrow line region (e.g Cameron et al 1993, Braatz et al 
1993, GaUiano et al. 2005). Efstathiou, Hough & Young 
(1995) presented a model for the nuclear spectrum of 
NGC1068 that involved a tapered disc (Efstathiou & 
Rowan- Robinson 1995) and optically thin dust in a con- 
ical region. The latter may actually reside in the narrow 
line region (NLR) clouds. The model considered in this 
letter for the AGN contribution to the spectrum of FSC 
10214+4724 is a combination of the tapered disc mod- 
els of Efstathiou & Rowan- Robinson (1995) and a sim- 
ple model for the emission of clouds that are associated 
with the NLR. A tapered disc has a thickness that in- 
creases linearly with distance from the central source in 
the inner part of the disc but assumes a constant height 
in the outer disc. Although the models of Efstathiou & 
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Rowan-Robinson do not consider clumping of the dust 
in clouds they can be considered as a reasonable approx- 
imation to the dust distribution in AGN if the distance 
between clouds is much less than the size of the torus. 



3 RESULTS AND DISCUSSION 

Detailed modeling of the emission of NLR clouds that 
are illuminated by a central source is beyond the scope 
of this letter. The emission of such clouds is approxi- 
mated with that of spherical clouds of dust at a con- 
stant temperature. It is assumed that the density in the 
cloud is uniform and the optical depth in the V band 
from the center to the surface is 10. The same grain 
mixture used for the dusty starburst radiative transfer 
calculations is assumed. A good fit to the Spitzer spec- 
trum can be obtained by assuming that the NLR dust 
is concentrated in clouds at a temperature of 200 and 
610K with the 200K cloud(s) being about an order of 
magnitude more luminous. This is broadly in agreement 
with the findings of Teplitz et al. (2006). The starburst 
model assumed is the tv ~ 200 model shown in Fig. 1. 

The model assumed for the torus is the most geo- 
metrically thin of the models considered by Efstathiou 
& Rowan- Robinson (Oi = 30°). This implies a torus 
half opening angle of 60° which, however, is very poorly 
constrained. The model also assumes a ratio of inner to 
outer disc radius of 0.05, a ratio of height to outer ra- 
dius of 0.5 and an equatorial optical depth at lOOOA of 
1000. It is assumed that the tapered dusty disc in FSC 
10214-1-4724 is viewed edge-on. This means that its mid- 
infrared radiation is heavily suppressed. The sum of the 
emission of the NLR dust and the disc is shown in Fig. 2 
and 3 whereas the individual spectra are shown in Fig. 
4. Note that because the starburst is very optically thick 
it does not make a significant contribution to the mid- 
infrared flux. This can explain why there is no evidence 
of PAH features in the Spitzer spectrum published by 
Teplitz et al. except for a weak 6.2/im feature that is 
reproduced by the model. 

As can be seen from the plot of the edge-on and 
face-on spectra emitted by the torus in Fig. 4, the AGN 
luminosity is highly anisotropic. To obtain the total lu- 
minosity one needs to multiply the 'apparent' torus lu- 
minosity with the anisotropy correction factor A which 
is defined to be 

r^^S{i') sini' di' 

where S{i) is the bolometric emission at angle i. For the 
particular tapered disc model and inclination, A is found 
to be 18.1. The emission by the NLR dust is assumed 
to be isotropic. With this correction the AGN luminos- 
ity is predicted to be 6.9 x W^^Lq. The luminosity of 
the 200K NLR dust is predicted to be 1.0 x IO^Lq or 
15% of the total AGN luminosity. As mentioned above 
the luminosity of the 610K dust is about an order of 
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Figure 2. Fit to the spectral energy distribution of FSC 
10214-1-4724. The dotted and dashed lines show the contri- 
butions of the AGN and starburst respectively to the total 
emission (solid line). Data from Teplitz et al (2006), IRAS, 
Rowan-Robinson et al. (1993). 

magnitude lower or 2% of the total. As the NLR clouds 
are assumed to be optically thick, the fraction of total 
AGN luminosity emitted by the 200 and 610K dust is 
also the covering factor of the NLR dust. The starburst 
luminosity is predicted to be 7.6 x W^^Lq. Note that 
all of the quoted luminosities need to be reduced by a 
factor 10-100 due to gravitational lensing. Lens models 
suggest that the optical emission may be more highly 
magrufied than the infrared emission that dominates the 
bolometric luminosity. Eisenhardt et al. (1996) estimate 
a magnification factor for the infrared emission of 30 
which implies that the AGN luminosity is of the order 
of 2.3 X 10^^ Lq. Note, however, that since the emission 
by the NLR dust and the torus does not arise from the 
same region it is possible that their magnification will 
be different. 

As it is demonstrated in Figs. 2, 3 and 4 the pres- 
ence of dust in the narrow-line region of an AGN nu- 
cleus, with a covering factor of 17%, can explain the 
fact that silicate features appear in emission both in lu- 
minous type 2 and type 1 AGN. The presence of these 
features in type 2 AGN is entirely consistent with the 
unified scheme (Antonucci 1993) as a dust torus that 
makes an insignificant contribution in the mid-infrared 
may still be present. The anisotropy of the emission from 
such a torus is so high that when the same system is 
viewed face-on the torus emission dominates and the 
effect of the presence of the NLR dust is to make the 
silicate emission features more prominent. The face-on 
spectrum shown in Fig. 4 is similar to the spectra of 
quasars recently observed with Spitzer. 

A question that naturally arises is why lower lumi- 
nosity type 2 AGN do not show the same effect. There 
are probably two reasons for this; firstly in lower lumi- 
nosity type 2s the ratio of the luminosity of the acom- 
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Figure 3. As Fig. 2 but siiowing tlie fit to tlie mid-infrared 
spectrum in more detail. 
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Figure 4. Comparison of tlie edge-on and face-on views of 
FSC 10214-1-4724 according to the model presented here. The 
thick (red) dotted line gives the spectrum emitted by the 
edge-on torus and the dashed line (magenta) gives the emis- 
sion of the dust in the cones. The thick (red) solid line gives 
the total emission in the edge-on case. The thin (blue) dot- 
ted line gives the predicted spectrum for a face-on view of 
the torus and the thin (blue) solid line gives the sum of that 
emission to the emission of the dust in the cones which is as- 
sumed to radiate isotropically. The vertical scale is arbitrary. 

panying starburst to that of the AGN is higher so what 
we see in the mid-infrared is the complex of PAH fea- 
tures and silicate absorption features emitted by the 
starburst. Secondly the torus opening angle in lower lu- 
minosity AGN may be smaller and therefore the frac- 
tion of the source luminosity that is intercepted by NLR 
dust is smaller. This was already suggested by Lawrence 
(1991) who estimated the statistics of type 1 and type 2 
AGN from different surveys. As proposed by Lawrence 



this dependence on luminosity may arise if the thickness 
of the torus is independent of luminosity and its inner 
radius is controlled by dust sublimation. 
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